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ABSTRACT 

We report on CO (J = 2 —>■ 1) mapping with the IRAM 30-m HERA receiver array 
of CGCG 97-079, an irregular galaxy in the merging galaxy cluster Abell 1367 {z = 
0.022). We find that ~ 80% of the detected CO (J = 2 —>■ 1) is projected within a 
16 arcsec^ (6.5 kpe^) region to the north and west of the optical/NlR centre, with 
the intensity maximum offset ~ lOarcsec (4kpc) NW of the optical/NIR centre and 
~ 7 arcsec (3 kpc) south-east of the Hi intensity maximum. Evolutionary synthesis 
models indicate CGCG 97-079 experienced a burst of star formation ^ 10® yr ago, 
most likely triggered by a tidal interaction with GGGG 97-073. For GGCG 97-079 we 
deduce an infall velocity to the cluster of ^ 1000 kms“^ and moderate ram pressure 
(Pram ~ 10“^^ dyncm“^). The observed offset in GGGG 97-079 of the highest density 
HI and GO(J = 2^1) from the stellar components has not previously been observed 
in galaxies currently undergoing ram pressure stripping, although previous detailed 
studies of gas morphology and kinematics during ram pressure stripping were restricted 
to significantly more massive galaxies with deeper gravitational potential wells. We 
conclude the observed cold gas density maxima offsets are most likely the result of 
ram pressure and/or the high-speed tidal interaction with GGGG 97-073. However 
ram pressure stripping is likely to be playing a major role in the perturbation of lower 
density gas. 

Key words: galaxies: CO — galaxies: ISM — galaxies:clusters:individual: (Abell 1367, 
CGCG 97-079) 


1 INTRODUCTION 


The consensus view is that, in low-redshift galaxy clus¬ 
ters, ram pressure stripping by the intra-cluster medium 
(ICM) is the dominant mechanism accelerating the evolution 
of late-type galaxies (Ivan Gorkoinl l2004l : iBoselli Sz Gavaz^ 
I 2 OI 4 I : IChung et al.l I 2 OO 9 I I. Models for cluster spirals sub¬ 
ject to ram pressure stripp i ng (e . g.. iRo ediger fc HensleJ 
AosI: Rogdigerfc Brii gg^ l2007l : iKapferer et al.l 20091 : 


Tonnesen fc Brvan I 2 OO 9 II indicate that removal of the inter- 


* E-mail: tom.scott@astro.up.pt (TCS) 


stellar medium (ISM) proceeds progressively during a spi¬ 
ral’s infall to the cluster. Initially only the loosely bound 
peripheral interstellar medium (ISM), principally Hi, is 
stripped. The higher column density Hi and molecular gas 
located in the deepest parts of a galaxy’s potential well, 
traced by the optical/Near Infared (NIR) intensity max¬ 
ima, remain unperturbed until late in the stripping pro¬ 
cess and requires high velocities relative to the ISM/ICM 
(Vrei) and high ICM densities encountered during a galaxy’s 
traverse of a cluster core in order to remove them (e.g., 
iRoediger fc Hensl^ I 2 OO 5 I : iRoediger fc Briigg^ 120071 ') . The 
truncated Ha and dust disks observed in Hi deficient spi- 
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dyne cm ^), as 
280 kpc from the 


rals are consistent with this picture of outside-in IShf re¬ 
moval (iKoopmann fc KennevI l2004l : ICortese et alj l2010h . 
Hi and CO observations and modelling of the Virgo spi - 
rals NGC 4522 and NGC 4330 bv IVollmer et al] ll2008l . l2012l l 
are examples showing that during ram pressure stripping 
the high column density Hi and molecular gas remain lo¬ 
cated deep within the gravitational well traced by the NIR 
and opticall y most luminous pa rt of the galaxies. Moreover 
NGG 484 8 Ivo llmer et al■ll2001^ in Coma {Lx =7.8 x 10^^ 
ergs s~^: IPIionis et al.l 20091 ) indicates that even in X-ray 
luminous clusters a spiral’s transit of the cluster core may 
remove Hi without causing a molecular gas def icienc:J^, al¬ 
though the recent work by iBoselli et al.l (l2014l i for a sam¬ 
ple of Virgo galaxies shows a c orrelation between m olecu- 
lar and Hi deficiencies (see also iFumagalli et al.ll200^ 1. The 
ISM content, morphology and kinematic signatures result¬ 
ing from extreme ram pressure 10 
is proposed for ESO 137-001 projected 
center of the Norma cluster {Mdyn ^ 1 x 10^^ Mfri) from mod¬ 
elling with Kei ~ 3000 kms“^ bv ijachvm et al.l ll2014li . are 
less well studied. Those authors argue extreme ram pres¬ 
sure stripping has already removed all detectable Hi from 
the disk of ESO 137-001 and is now in the process of strip¬ 
ping the molecular gas, as well as being partially responsi¬ 
ble for a molecular tail containin g approximately half of the 
galaxy’s observed molecular gas. Jachvm et al.l ll2014h note 
the stripped Hi is likely to have have been converted to 
other phases, including the observed molecular and X-ray 
tail emission. Even ESO 137-001 appears to have followed 
the sequential removal of increasingly dense ISM from the 
galactic disk during ram pressure stripping, with stripping 
of the molecular disk being delayed until the final stages 
which are now being witnessed. 

There have been relatively few detailed studies of 
changes in H l and CO morphology and gas kinematics dur¬ 
ing ram pressure stripping in gas-rich late-type spirals and 
in particular of dwarf irregular galaxies. Although a Virgo 
Irr dwarf, IC 3418, with about twice the stellar mass of the 
SMC, contains evidence of having lost almost all its Hi and 
molecular gas by ram pressur e stripping within the last few 
X 10® yr jKennev et al.ll2014l l. A sample of late-type dwarf 
galaxies near the centres of clusters, which are presumed 
to have lost their Hi by ram pressure stripping, display 
truncated Ha disks (iFossati et al.l 12013 1. Also we see clear 
evidence of truncated UV star-forming disk s in both spi¬ 
ral a nd dwarf galaxies in the Virgo cluster (|Cortese et al.l 
I 2 OI 2 II . These studies imply dwarfs are subject to the same 
sequentia l outside-in gas stripping as higher mass late-type 
galaxies JBoselli fc Gavaz^ l2014^ . a proposition which is 
supported by the detection of CO in the disk of IC 3418. 

The scenario depicted above is challenged by our recent 
studies of late-type galaxies in the spira l-rich merging clus - 
ter Abell 1367 (Lx = 1-25 x 10"^"^ e gs s~^ :l Plionis et al.l 200^1 
which has a Mdyn ~ 6.9 x 10^"^ Mq (iBoselli fc Gavazzi2006li . 
about half that of Coma or Norma, and hosts two sub¬ 
clusters in the early stages of an a pproximatel y equa l mass 
merger (iDonnellv et al.l Il998l i . In Iscott et al.l (120101 . here- 


^ A spiral’s H 2 or HI deficiency is defined as the log of the ratio 
of the expected to observed gas mass. Negative values indicate an 
excess. 


after Paper I), we found an unusually high frequency of off¬ 
sets between the H l and optical intensity maxima in a sam¬ 
ple of the cluster’s late-type galaxies. Our multi-pointing 
CO (J = 1 —t 0) and CO (J = 2 —» 1 ) observations with the 
IRAA0 30-m telescope dScott et al.ll20l3 . hereafter Paper 
H), provided indications of perturbed molecular disks and 
molecular gas excesses in a subset of the Paper I sample. 
Together, these papers raise the question whether a mecha¬ 
nism other than standard ram pressure stripping is driving 
the evolution of spiral galaxies in A 1367, and in merging 
clusters in general. 

CGCG 97-079 is an optically-irregular galaxy with D 25 
= 45 arcsec (18.6 kpc), an estimated M* ~ 1.2 x 10® Mq, 
(i.e., of order of the Large Magellanic Cloud or M33) and is 
projected 420 kpc to the NW of the NW sub-cluster core of 
A 1367 (Figure 1 in Paper II). That figure also shows the X- 
ray emission (ROSAT) from the A 1367 ICM. Its 50-100 kpc 
ionised gas and radio continuum tails oriented away from the 
cluster centre (Figure [T]), provide evidence that CGCG 97- 
079 is undergoing ram pressure stripping (ICavazzi fc Jaffel 
1 19871 : [Cavazzi et al.ll200l] : IBoselli fc Gavazziil2014 '). However 
our previous observations (Papers I and H) suggest its neu¬ 
tral (atomic and molecular) gas distributions are inconsis¬ 
tent with current ram pressure stripping models. In partic¬ 
ular its molecular and atomic gas intensity maxima lie ~6 
kpc and ~12kpc, respectively, to the NW from its optical 
intensity maximum as derived from sparsely sampled CO 
maps and low resolution H l data. CGCG 97-079 is the clear¬ 
est case in A 1367 of Hi and CO intensity maxima offsets, 
in apparent contradiction to ram pressure models, making 
it an interesting target for the investigation of mechanisms 
driving the the evolution of ISM in galaxies in A 1367. 

This paper reports on follow-up CO (J = 2 ^ 1) map¬ 
ping of CGCG 97-079 with the IRAM 30-m HERA receiver 
array. We also utilised Hi C-array data from the NRAC0 
VLA archive (project AG264), a Spitzer IRAcQ 4.5 pm im¬ 
age (ID: AORKEY25789696 taken from the Infrared Science 
Archive IRSA), a GALEX F\JY image (A1367.SPECfrV-fd- 
int.fits) retrieved from the MAST archive 0, and an Ha im¬ 
age from GOLDMimEI 

Section [ 2 ] gives details of the HERA observations, with 
observational results in section [S] A discussion follows in 
section |4] with concluding remarks in section [5] Based on 
a redshift to A1367 of 0.022 and as suming = 0.3 , 
Ha = 0.7, and Ho = 72kms~^ Mpc“^ (ISoergel et akllioOTh 
the distance to the cluster is 92 Mpc and the angular scale 
is 1 arcmin ~ 24.8 kpc. All equatorial positions referred to 
throughout this paper are in J2000.0. References to the opti¬ 
cal centre are to the position ll''43™13!4, -|-20°00'17" from 
the NASA/IPAC Extragalactic Database (NED). 


^ Institut de radioastronomie millimetrique (IRAM) is supported 
by CNRS/INSU (France), the MPG (Germany) and the IGN 
(Spain). 

® The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 

InfraRed Array Gamera. 

® Mikulski Archive for Space Telescopes. _ 

® Galaxy Online Database Milano Network. llGavazzi et al.ll2003l . 
http://goldmine.mib.infn.it/). 
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U‘'43”15= 43™I0® 43“05* 43”00* 42“55= 

Right Ascension (J2000) 


Figure 1. CGCG 97-079: SDSS r-band image with the black (shorter ) and white (longer) arrows indicating, respectively, the approximate 
extent and orientation of the radio continuum llGavazzi &: JafF3ll987li and Ha tails of CGCG 97-073 and CGCG 97-079. The white arrow 
head of the CGCG 97-079 tail marks the approximate position where the Ha tails from deep imaging from Boselli &; Gavazzi (2014) 
appear to meet. Ha contours (white) are from a GOLDMine image. The position of the nearest neighbouring galaxy, A 1367 GP82 1236 
is also marked. 


2 OBSERVATIONS-HERA CO MAPPING 

(J = 2 —I 1) emission line mapping of CGCG 97-079 
was carried out with the HERA 3x3 multi-beam receiver 
array on the IRAM 30-m telescope at Pico Veleta, Spain, in 
slightly under sampled dual polarisation “On the Fly” mode 
with position switching. Receivers used were HERA 1 and 2 
with WILMA 1 and 2 backends (1024 channels with 2 MHz 
separation). Further observational parameters are given in 
Table [T] The total integration time was 13 hr which was 
accumulated in March 2010 and January 2012. CO (J = 2 —4 
1 ) was mapped ina~ 1.1x1.1 arcmin^ region of sky centred 
at (-10.6, -1-10.6 arcsec) from the galsixy’s optical centre. 
After every three scans (~ 15 minutes) a chopper wheel 
calibration was carried out. Pointing was checked at between 
1 and 2 hour intervals using the broadband continuum of 
3C273. 

Data reduction was carried out using GILDAS software 
and excluded spectra taken under particularly poor condi¬ 
tions. Spectra were then summed position by position. The 
data cube was produced using the PLAIT task and excluded 
unevenly-sampled edge regions. Observations covered the 
velocity range 6552 kms“^ to 7425 kms“^ with a velocity 
resolution of 11 kms“^. The data cube was blanked using 
the AIP^ software package task blank resulting in a final 
cube containing 42 channels binned to a velocity resolution 
of 21kms“^ with an angular resolution of 11.5 arcsec. 


Astronomical Image Processing System 


Table 1. HERA observational parameters 


Rest frequency 

(J = 

2 l)[GHz] 

230.5 

Sky frequency 


[GHz] 

225.3 

Primary Beam 

FWHP_^ 

[arcsec] 

11.5 

Primary Beam 

FWHP 

[kpc] 

4.5 

Receivers 

HERA 


1 & 2 

Backend channel width 


[MHz] 

2 

Tsys 

T @ 225 GHz 


[K] 

300±100 

0.3±0.1 

Expected pointing accuracy 

Tmb 

[arcsec] 

2 

1.73 Ta* 


“ Full width half power 

3 OBSERVATIONAL RESULTS 

The GO (J = 2 —>■ 1) contours from the HERA integrated 
intensity map, in Figure [2] (top left), show that the bulk of 
CO (J = 2 ^ 1) was detected within a continuous emission 
region approximately 16 arcsec^ (7kpc^) in area, projected 
west and north of the optical centre. This continuous emis¬ 
sion region is referred to as the “CO disk”. The CO disk 
contains the CO (J = 2 —^ 1) intensity maximum (CO max¬ 
imum), indicated with yellow plus signs in Figure [2] and is 
located 9.9 arcsec (4.1 kpc) NW of the Spitzer 4.5 /rm 
intensity maximum {PA = 287°). Also marked in the same 
figure is a detached CO (J = 2 —> 1) clump (Cl) ~ 29 arcsec 
(12 kpc) NW of the optical centre. Cl is detected in 3 con¬ 
tiguous channels at > 3 ct. We are cautious about declaring 
Cl a detection because of its low peak signal to noise (S/N) 
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of 3.9, its lack of a counterpart at any other wavelength 
and its proximity to the map edge. However on balance its 
CO (J = 2 —>■ 1) spectrum (Figure [3]- bottom left), CO 
(J = 2 —>■ 1) single beam detection from Paper II, velocity 
and the total H 2 mass including Cl (see next paragraph) 
support Cl being a real feature. The CO disk morphology of 
CGCG 97-079 differs dramatically from t hat of local spirals 
dNishivama et al]l200ll : lLerov et al.ll2009h . If it followed that 

of the HERACLES nearby spiral sample ijLerov et al.ll2009l l. 

it would display an exponential distribution with its origin 
at the optical centre with a typical expected scale length of 
0.23r 25 (i.e., 3.4arcsec or 1.4kpc). 

L(CO) 2 -i from the CO disk and Cl is 1.4±0.1 x 
10® K kms“^pc^ (Tmb scale) which converted to an H 2 
mass, using equation Bl][f| from Paper H (X factor of 
2.62 X 10^'^ molecules cm~^/Kkms“^), gives M{H 2 ) = 0.8 
±0.1 X 10^ Mq. This compares to the earlier H 2 mass es¬ 
timate based on single beam CO (J = 1 —>■ 0) pointings 
(Paper H) of 1.07 xlO® M©. Figure [2] - top left shows 
the position of the CO (J = 1 ^ 0) full width half power 
(FWHP) single beam pointings (dashed yellow circles), from 
which the earlier H 2 mass was estimated. The HERA data 
presented here have a higher signal to noise. A comparison 
of spectra extracted from the current data with the single 
CO (J = 2 —^ 1) pointings from Paper H leads us to suspect 
the Paper II masses are less reliable. Further properties of 
the CO disk and Cl are given in Table[2]with Vhei, and W 20 
calculated using the methods from Paper H. 

Figure [ 3 ] (top left) shows the CO (J = 2 —>■ 1) velocity 
field for the CO disk and Cl in the regions where emission 
is > 3 (j in the integrated intensity map. From this velocity 
field and PV diagram (Figure |4] - top left) we see a sys¬ 
tematic increase in velocity from 6960 kms“^ along an axis 
joining the optical centre and the CO (J = 2 —>■ 1) inten¬ 
sity maximum where the velocity reaches 7020 kms“^. From 
the CO (J = 2 —>■ 1) maximum the velocity remains more 
or less constant between 7020 kms“^ to 7040 kms“^ out to 
the Cl clump. A second velocity gradient is seen in the PV 
diagram, Figure[4](top right), running NE (6936 kms“^) to 
SW (7085 kms“^) along an axis approximately aligned with 
knots D and knot E {PA ^ 244°). Knot i dentif ication letters 
follow the Hq knots from lCavazzi et ^ lll995ll as marked in 
Figure [5] (bottom left). This second PV diagram shows two 
velocity components, one at 7000 kms“^ and the other 
at 7050 kms“^. 

The HERA CO (J = 2 —^ 1) spectra for the CO disk, 
Cl and a combined CO disk ± Cl spectrum are shown in 
Figure O The peak signal to noise (S/N) in the spectra are 
6.8, 3.9, and 7.3, respectively. The CO (J = 2 —^ 1) velocity 
of the CO disk ± Cl (V = 7030 ±12kms“^) agrees within 
the errors with the Hi velocity of 7019 ±21 kms“^ (Paper 
I) as does the HERA CO (J = 2 -i- 1) W 20 (240±25 kms“^) 
and Hi W 20 (216 ±41 kms“^), where W 20 is the full width 
of the line at 20% of the peak emission in the spectrum. 


® Equation Bll requires L(CO)i_o and to apply the equation 
we used the conversion: R21/R10 = 0.75 


4 DISCUSSION 


In the following we summarise the properties of the gas 
and stellar components of CGCG 97-079, based on available 
data, incl uding a series of papers by Gavazzi and collabora- 


tors (e.g.. Gavazzil 19781: 

Gavazzi & Jaffd 19871: Boselli et al.l 

ll994lGavazzi et al.lll995 

.l200ll:lBoselli & Gavazzill2014lL We 


subsequently consider the scenarios most likely to explain 
the accumulated observations in the concluding remarks 
(section [S]) 


4.1 The old stellar component 


A combination of perturbed ISM and an unperturbed old 
stellar population is the key signatu re of an ongoing ram 
pressure stripping interaction (e.g., iKennev et al.l l2004l L 
Hubble Space Telescope imaging of E SQ 137-001 which is 
suffering strong ram pressure stripping lljachvm et al.|[^014h 
shows narrow trails of blue optical emission, aligned down¬ 
stream of the optical disk, as well as extra-planar dust 
above the inner part of the optical disk. So it is clear that 
extra-planar emission/absorption from dust mixed with the 
stripped gas as well as emission from young stars can result 
from ram pressure stripping. In general infrared (IR) emis¬ 
sion is a better tracer of the old stellar population than opti¬ 
cal emission. But IR emission from a galaxy’s old stellar pop¬ 
ulation can be contaminated by dust re-emission following 
absorption of UV photons from recently formed high mass 
stars and stellar emission fro m RGB stars with ages of 10 to 
100 Mvr. lBoselli et all (l2004h give the dust emission contam¬ 
ination at 6.75 /rm as up to 8 0% for Sc gal a xies a nd 50% for 
blue compact dwarfs. But as iBoselli et al.l (|2004l 'l note MIR 
dust re-emission is more intense in massive quiescent galax¬ 
ies than in star forming metal poor ones. CGCG 97-079 is 
both a dwarf and metal poor {FejH — -0.62 lMouhcine et al.l 
[ 2 ^. However, while contamination by dust re-emission is 
significant at mid infrared (MIR) wavelengths (including the 
6.75 and 15 /rm ) it is minimal in th e Syitzer 3.6/xm an d 
4.5/rm llFazid I2O O5I: lUtomo et alJ I20l4 iMeidt et al'l 120141 ') . 
Modelling bv lUtomo et al.l 'i 2014 1 indicates the contamina¬ 
tion (including from polycyclic aromatic hydrocarbon) re¬ 
mains minimal (of the order of a few %) at 3.6/im and 
4.5/rm even for galaxies with Ha equivalent widths {EW) 
similar to CGGG97-079 (EIV(Ha) =129 per GOLDMine). 
For a disk galaxy contamination from young stellar emission 
(principally from red supergiants) at ~2.3/im is estimated 
at ~3%, but in star forming regions w ithin the galax y the 
contamination can be as high as 33% (|Rhoadslll998^ . The 
strength of the emission from red supergiants in the wave¬ 
length range of 2.3 /rm and 4.5 /rm declines with increasing 
wavelength llBaron et al.ll20f3 ). This suggests that red su¬ 
pergiant contamination at 4.5 /rm is likely to be less than at 
2.3 /rm by a factor of a few. Because of the non-availability 
of a Spitzer 3.6pm image we were unable to make the poly¬ 
cyclic aromatic hydrocarbon (PAH) and young stellar emis¬ 
sion contam ination corrections developed by llMeidt et al.l 
|2012| . I 2 OI 4 I I . We conclude that in the absence of these cor¬ 
rections the 4.5pm band image provides the best available, 
tracer of the old stellar pop ulation. This ban d is known to 
suffer some CO absorption (iMeidt et al.ll20l3 L 

Figure[2] (bottom right) shows the Spitzer 4.5 pm image 
for CGCG 97-079, with contours from the same image with 
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Table 2. CO (J = 2 —> 1) properties from HERA 


Property 

Unit 

CO disk 

Cl 

Disk ± Cl 

L(CO)2-i [T,„b] 

[K kms“^pc^] 

1.2 ±0.1 X 10“ 

0.2 ±0.05 X 10“ 

1.4 ±0.1 X 10^ 

M(H2) 

[10® M©] 

0.7 ±0.1 

0.1 ±0.03 

0.8 ±0.1 

Vhel 

[kms“^] 

7000 ±13 

7045 ±31 

7030 ±12 

W 20 

[kms“^] 

220 ±26 

210 ±62 

240 ±25 

rms per channel [Tmb] 

[mK] 



3.1 

rms integrated map [Tmb] 

[K kms-l] 



0.23 


a 3-pixel wide boxcar smoothing applied. The un-smoothed 
image had a pixel size =1.2 arcsec and resolution ~ 2.5 ar- 
sec). Two principal intensity maxima are seen in the figure: 
an infrared counterpart to Ha knot A, and a more elongated 
counterpart approximately coinciding with Ha knots B and 
C. The 4.5 /rni intensity maximum is projected at approxi¬ 
mately the same position as the optical centre and knot B. 
The figure shows the intensity of the 4.5 /rm emission falls 
steadily south of the optical centre, but to the N and NW of 
it there is an extensive region of diffuse emission. An uncata¬ 
logued 4.5 /rm feature marked F on the figure is not present 
in the 2MASS, optical (SDSS) or Ha images, although it 
does have a counterpart in our J - band image obtained 
with the SPA0 2.1 m telescope. Knot F’s colour suggests it 
is a higher redshift background galaxy. The optical counter¬ 
parts of the two objects in the same figure ~ 10 arcsec N 
of knot F are catalogued as 23rd magnitude stars in SDSS. 

Depending on the weighting given to the high and low in¬ 
tensity 4.5 /rm emission there are several possible alternative 
inclinations and morphologies for the old stellar disk. Three 
possible old stellar disk orientations are over plotted as el¬ 
lipses in Figure O The yellow ellipse represents an edge-on, 
~ 90° inclination disk {PA = 275°) aligned with knots A, 
B, C and D, the dashed white ellipse is a lower inclination 
(60°) disk {PA = 286°) aligned with the outer edges of the 
detected 4.5 /rm, and the solid white ellipse represents an 
edge on, ~ 90° inclination disk {PA = 286°) aligned with 
knots B and C. In case of the high inclination disk options 
the diffuse 4.5 /rm emission to the north of the optical cen¬ 
tre, and particularly the upward curve at the western disk 
edge, would be evidence of a tidal perturbation. On the other 
hand for the lower inclination (60°) disk the northern dif¬ 
fuse emission would just be part of the disk. Because of these 
and further alternative interpretations of the 4.5 fim data, 
as well as our inability to make the PAH and young stel¬ 
lar emission contamination corrections, we were unable to 
definitively determine whether the morphology of the old 
stellar disk is perturbed or not. 


4.2 HI and CO gas components 


CGCG97-079 is particularly gas rich, with M(H 2 -I- Hi -|- 
He) = 2.9 X 10® M© (Mm = 1.3 x 10® M© - P^er I) 
and taking M* = 1.2 x 10® M© gives a gas fractioiT®!. /gas 
= 0.7. Its M(H l)/M« = 1.25 is consistent with that found 
for dwar f (IZhan g et al.ll2012| j and smaller late-type galaxies 
dCortese et al. 2011 1 like the LMC and M33. Based on the 


® San Pedro Martir, Mexico (Observatorio Astronomico Na- 
cional) 

Ratio of gas mass to the sum of gas and stellar masses. 


HERA observations the H 2 excess (including Cl) is -0.50, in 
good agreement with the H 2 excess from Paper H (-0.64). As 
discussed in Paper H, there are significant uncertainties both 
in the determination of the H 2 mass of individual galaxies 
and the calibration of H 2 deficiencies from samples of unre¬ 
solved CO observations, but based on the available data we 
can, at a minimum, co nclude that CGCG 97 -079 is not de¬ 
ficient in molecular gas. ICavazzi et al.l (1200 ill estimated the 
mass of Ha in the ionised tail as 9.6 x 10® M©, which could 
account for almost all of the Hi deficiency (0.25) reported 
for CGCG 97-079 in Paper I, provided the Ha emission is 
from gas that was originally Hi stripped from the galaxy. 
The NW offset of the CO (J = 2 —^ 1) distribution from the 
optical centre, presented in section [3l is further illustrated 
in Figure |4] (lower panel), which shows the normalised CO 
(J = 2 —>■ 1) emission integrated along the axis of knots B 
and C in comparison with 4.5 /rm emission along the same 
axis. The displacement of the CO (J = 2 —>■ 1) from the 
highest density old stellar component of the galaxy is clearly 
seen in the figure. 

Figure [2] (top right) shows that the H l intensity max¬ 
imum from the VLA C-array is offset by 7 arcsec (3 kpc) 
north-west of the CO (J = 2 —>■ 1) maximum and has a 
column density derived from the VLA data (C-|-D-array) 
of 6.6 X 10®® atomcm“® (iHota fc Saikia|[2007l f . This offset is 
consistent with our deeper VLA - D array observation in 
Paper I. The Hota VLA data (C-l-D-array) map (their fig¬ 
ure 14) shows the Hi velocities increasing from 6923 kms“^ 
in the NE to 7117 kms“^ in the SW along the knot E to 
D axis. A similar CO (J = 2 ^ 1) velocity gradient is ob¬ 
served along the same axis and the CO (J = 2 —>■ 1) PV 
diagram for a cut along this axis indicates there are two ve¬ 
locity components (Section |3]). Also detailed in section [3] is 
the CO (J = 2 —>■ 1) velocity gradient between the optical 
centre and the CO (J = 2 —> 1) maximum with a decrease 
in the gradient NW of the CO (J = 2 —>■ 1) maximum. 

Both the offsets of the H l and CO (J = 2 —> 1) max¬ 
ima from the optical/4.5 /im centre and the segregation of 
the high column density Hi and CO (J = 2—^1) maxima 
contrast with NGC4522 and NGC4330 in the Virgo cluster 
dVollmer et al.|[2008l . 1201^ 1. where the high column density 
Hi and CO remain projected at the optical centre during 
ram pressure stripping. Although NGC4522 and NGC4330 
are both significantly more massive than CGCG 97-079. 


4.3 Recent star formation 

Ha (Figure [2] bottom left) and FUV GALEX (Figure [6]l 
emission trace star formation on time scales of 10^ yr and 
10® yr respectively (iBoselli et al.l l2009l j and are both con¬ 
centrated east of the CO (J = 2 —> 1) maximum. About 
30% of the CGCG 97-079 FUV emission is detected within 
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Figure 2. CGCG 97-079 top left: CO (J = 2 —>■ 1) HERA integrated intensity map contours (gray) overlaid on a composite SDSS false 
colour g,r,i - band image. The first contour level is at 3 cr, equivalent to 0.68 K kms“^ (Tnib) and subsequent contours are at 1.0, 1.3, 
1.7, and 2.0 K kms“^ (Tmb)- The white circle indicates the size of the CO (J = 2 —> 1) beam and Cl indicates the position of the 
detached clump. The IRAM 30-m 22 arcsec FWHP single beams for the CO (J = 1 —> 0) Paper I observations at the optical centre 
and offset RA = -15.6 arcsec, DEC =15.6 arcsec are indicated with large dashed yellow circles. Top right: Hi (VLA C-array) contours 
(cyan) overlaid on a 4.5 /rm image (resolution ~ 2 arcsec). The contours are at 2, 2.5, 3, and 3.5 cr. The white circle indicates the size of 
the VLA C-array beam. Bottom left: Ho contours (green) from GOLDMine (resolution ~ 3 arcsec) overlaid on a composite SDSS false 
colour g,r,i - band image. The contours are at 20, 50, 100, 200, 400 and 800 cr. The letters A to E mark the 5 principal HII regions (Ho 
knots), following the naming convention from Figure [2] of llGavazzi et al.iri995h . Bottom right: Spitzer contours from a 4.5 /cm boxcar 
smoothed image, overlaid on a 4.5 /cm un-smoothed image. The contours are at 4.5, 5.5, 10, 15, 25, 35, 45, 55, 65 and 75 cr with the 
4.5 cr contour equivalent to 0.085 MJy sr“^. The resolution of the 4.5 fim image is ~ 2.5 arcsec. The orange asterisk marks the optical 
centre and the cyan ’X’ in box marking the position of the HI intensity maximum from the VLA C—array. The red ’X’ corresponds to 
the 4.5 /tm maximum and the yellow plus symbol the GO (J = 2 —>■ 1) maximum from the HERA map. 


a radius of 3 arcsec (1.25 kpc) of knot A. The Ha emis¬ 
sion is concentrated in knots A, B and C. These three knots 
are also the source of the st rongest 1.5 GHz ra dio contin¬ 
uum emission (Figure lb in iGavazzi et al.lll995l) . Figure [7] 
shows a. g — i band image produced from SDSS g and i 


band background subtracted images overlaid with HERA 
CO (J = 2 —1 1) contours. From the figure we see an ex¬ 
tensive area between knots B and E with g — i ~ 0.7. The 
figure shows this area correlates well with the high column 
density CO (J = 2 —l 1) emission. This strongly suggests 
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Figure 3. Top left panel: CO {J = 2 —>■ 1) velocity map with contours (kms“^) overlaid on a composite SDSS - g,r^i band image. The 
white circle at the bottom left of the panel shows the approximate size of the HERA CO (J = 2 —^ 1) beam. The colours and symbols 
are as in Figure [2] except the yellow plus symbol at the top right which marks Cl. Top right panel: HERA CO (J = 2 —> 1) spectra for 
the CO disk + Cl. Bottom left panel: HERA CO (J = 2 —>■ 1) spectrum for the CO disk. Bottom right panel: HERA CO (J = 2 —>■ 1) 
spectrum for Cl. The T^b for each channel in a spectrum (21 kms“^ velocity width per channel) is the sum of the Ta* for all pixels 
in the channel, converted to the T^b scale. The horizontal line at the bottom of each spectrum indicates the spectrum’s W 20 velocity 
width and the vertical cross bar the Yhel- 


the g — i value in this region is the result of dust extinction 
and the bluer g — i colours for knots A and D are partially 
the result of lower extinction compared to knots B and C. 
Even allowing for this, knot D appears to be dominated by 
a younger population of stars than the four other knots. An 
interesting feature in the Ha image is the ridge of Ha emis¬ 
sion along the knot E to D axis {PA = 244°) which has less 
well defined optical (SDSS), 4.5 ^m, H l and CO (J = 2 —>■ 1) 
counterparts. 


4.4 High impact low velocity tidal interactions 
with nearby neighbours 

Here we consider whether there is a neighbouring galaxy or 
galaxies which could have tidally disrupted CGCG 97-079 
within the last few x 10® yr assuming a typical group ve¬ 


locity of 250 kms“^. To try to identify any such galaxy 
we used a modified version of the tida l force (Q) and num - 
ber density {g — 2.5) parameters from IVerlev et aP (l2007h . 
applied to data from SDSS DR8. A Q < 0 indicates the 
gravitational forces within the galaxy exceed the gravita¬ 
tional force exerted by the neighbours. We find Q = —1.10 
for the 27 n eighbours in the wes tern sub-cloud of the NW 
sub-cluster (ICortese et al.l [20041 ) and Q — —1.25 for the 
nearest neighbour, A 1367 GP82 1236 (see Figure [T] at RA = 
11*'43“ 10.961% and DEC = -b20°0T47.0l"). The tidal per¬ 
turbation parameteJ*^ pgg — 0.004 for CGCG 97-079 with 

P 99 = Mgai and Mcomp are the masses 

of the galaxy and compani on respectively, d is the separation and 
r is the galaxy disk radius ijBvrd Valton^ll990lh where values 
of Pgg >0.1 likely lead to tidally induced star formation. 
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Figure 4. Top left panel: PV diagram (PA = 286°) centred at the optical centre. The angular offset, in arcsec, is negative to the SE 
(left) and positive to the NW (right). Top right panel: PV diagram (PA =244°) along the line connecting knot E and D centred at knot 
E. The angular offset, in arcsec, is negative to the NE (left) and positive to the SW (right). The horizontal arrows indicated the size 
of the HERA beam. Bottom left panel: CGCG 97-079 integrated Spitzer IRAC 4.5 pm (dashed line) and GO (J = 2 —>■ 1) (solid line) 
emission in normalised units along the PA = 286° centred on the optical centre of the galaxy. The horizontal axis shows the offset from 
the optical centre in kpc, with positive values to the NW. Bottom right panel: SDSS g - band image showing the orientation and lengths 
of the slices along which the PV diagrams where derived. 


A 1367 GP82 1236 is more than an order of magnitude be¬ 
low that expected to produce detectable molecular gas per¬ 
turbations in CGCG 97-079. We conclude from both the Q 
and pgg parameters and evidence that the perturbation oc¬ 
curred within ~ 10® yr (section 14.511 , that the substantial 
displacement of the high column density molecular gas in 
CGCG 97-079 is unlikely to be due to a low velocity high 
impact tidal interaction with A 1367 GP82 1236. 


4.5 High velocity interactions 

In this section we consider whether CGCG 97-079 may have 
undergone a recent high velocity (~ 1000 hms”*^) hydro- 
dynamic and/or tidal interaction. If the perturbation that 


caused H l and CO maxima offsets also drove the recent evo¬ 
lution of the stellar population we can infer the elapsed time 
since the interaction from the galaxy’s i ntegrated properties . 
GALEV evolutionary synthesis models jKotulIa et al.ll2009l) 
with realistic parameters produce a set of solutions for the 
observed M* = 1.2 x 10® Mq and SFR =1.06 Mq yr“^ (Pa¬ 
per II) for star bursts with burst strengths of ~ 0.20 - 0.25. 
The model constraint on these solutions from the optical 
and NIR colours is rather poor, using either the GALEV 
or FUV based Milky Way -1- internal extinction, but favour 
solutions with a time since the star burst of several x 10^ yr 
to a few X10® yr. The expected rela xation time scales for 
perturbed CO (J = 2 ^ 1) and Hi llHolwerda et al.ll201ll) 
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Figure 5. Stretched 4.5 fim band image of CGCG 97-079. The ellipses indicate three possible pre-interaction stellar disks: the yellow 
ellipse represents a edge on, ~ 90° inclination, disk (PA = 275°) encompassing knots A, B, G and D; the white dashed ellipse is a lower 
inclination (~ 53°) inclination disk (PA = 290°) tracing the outer edges of 4.5 pm — band emission; and the solid white ellipse represents 
an edge on, ~ 90° inclination, disk (PA = 286°) aligned with knots B and G. The cyan X in box marks the position of the Hi intensity 
maximum from the VLA G—array and the green plus symbols the GO (J = 2 —> 1) maximum and GI from the HERA map. The knot 
names are the same as used in Figure 2. 


together with GALEV modelling indicate the perturbation 
of the galaxy probably occurred ~ 10® yr ago. 

Th e deep Hq image of C GCG 97-079 and CGCG 97- 
073 in l|Boselli fc GavazMl2014l . their figure 21) shows each 
galaxy to have a spectacular Ha tail extending 100 kpc 
and 70 kpc from their respective optical disks and meet¬ 
ing in projection at ~ 11:42:57.7, -1-20:02:29.42 (see Figure 
[!}. At the vertex of the Ha tails the complex Hq mor¬ 
phology strongly suggests a high velocity tidal encounter 
occurred between the two g alaxies earlier i n the ir infall 
to the cluster as proposed in iGavazzi et ^ (1200 ll ). There 
are i ndications inlBo selli et ahn 19941) fat low significance) 
and ISivanandam et al.l 1201 'i) that the molecular gas in 
CGCG 97-073 is also offset in the direction of its tail. As¬ 
suming the most recent star burst in CGCG 97-079 was 
driven by the interaction at the tail intersection between 
7 X10^ yr to 2 x 10® yr ago (as indicated by the GALEV 
models), the CGCG 97-079 projected Ha tail length (~ 
4.3arcmin - 106 kpc), implies a velocity in the plane of the 
sky in the range of 500 kms“^ to 1500 kms“^ in agree- 
ment with the 1200 kms~^ velocity and 10^ ® yr derived by 
IGavazzi et al. Ilioil) for the survival time of the ionised tail. 

Simple modelling of ram pressure at CGCG 97-079’s 
distance from the NW sub-cluster core (Paper I) implies 
Pram « 10“^^ dyucm“^ for a Kei of 1000 kms“^. CGCG 97- 
079 remains gas rich (/gas = 0.7) with no overall gas defi¬ 
ciency (although marginally H l deficient) despite evidence 


from the Ha tails that ram pressure stripping has operated 
for least ~ 1 x 10® yr. This, the indication from section 
14. 2 1 that most of the stripped H l is now observed in the Ha 
tails and the estimate of Pram ~ 10“^^dyncm“^ all indi¬ 
cate CGCG 97-079 has suffered only moderate ram pressure 
stripping. Further support for moderate ram pressure comes 
from the absence of an H l counterpart to the Ha tail. Based 
on our VLA - D observations from Paper I we estimate there 
could at most be 3 x 10® M© of H l undetected in the tail, 
i.e. the Ha + potential undetected Hi in the tail account 
for the H i deficiency. 

Ram pressure stripping models have usually been 
applied to more massive and earlier Hubble-type spi¬ 
rals with classical bulges and lower gas fractions (e.g., 
spirals with M« = 3.8 x 10^° Mq and /gas = 0.08; 
Ifioediger fc Hensl^l2005l ). At Pram = 10“^^dyncm“^ the 
presence /absence of a bulge does not impact the stripping ef¬ 
ficiency JSteinhauser et aHl2012l ). But modelling of the ram 
pr essure stripping o f gas rich dwarves (M* = 0.6 x 10® M©) 
bv lMarcoIini et ^ ll2003l) have complete gas stripping time 
scales of between 2 — 4 x 10® yr for ram pressures at least an 
order of magnitude below that predicted for CGCG 97-079 
at Vrei = 1000 km s”*^. 

Confirmation of the ability of ram pressure stripping to 
remove almost all of the gas from a dwarf on timescales of a 
few X 10® yr comes from the study of IC 3418 llKennev et al.l 
l2014h . IC 3418 is a highly H l deficient Irr dwarf in Virgo pro- 
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Figure 6. CGCG 97-0979: Contours from a FUV (GALEX) im¬ 
age, overlaid on a composite SDSS false colour g,r,i - band im¬ 
age. The contours are at 4 ,15, 30, 45, 60, 100, 200, 300 and 400 
a with the 4 a equivalent to 1.74 xl0“^® erg sec“^ cm~^ A“^. 
The resolution of the FUV image is ~ 4 arcsec. 


jected 225 kpc from M87, with ~ l/3rd of the stellar mass 
of CGCG 97-079. The galaxy is thought to have undergone 
rapid ram pressure stripping within the last few x 10® yr 
which truncated its star formation. CO was marginally de¬ 
tected in the disk, which together evidence of super giant 
starformation 10® yr ago in the stellar disk, supports the 
proposition that the highest density gas remains within a 
dwarf’s disk until ram pressure stripping is complete. H l in 
IC 3418 was only detected within the starforming tail. 

Given CGCG 97-079 has a stellar mass, a gas fraction 
and other properties which are similar to dwarves as well as 
the evidence ram pressure stripping has operated for at least 
10® yr, we conclude that had CGCG 97-079 been subject to 
continuous ram pressure much stronger than predicted from 
Paper I (say Pram = ~ 1 x 10“^° dyncm“^) or its lower mass 
compared to generalised ram pressure models produced sig¬ 
nificantly more efficient Hi stripping, a substantial overall 
gas deficiency and H l tail would be observable by now. Since 
neither this deficiency nor an H l tail is observed, steady ram 
pressure seems unlikely to be the principal cause of the off¬ 
sets between the H l, CO (J = 2 —>■ 1) and stellar intensity 
maxima. However a sudden increase in ram pressure from 
a rapid inc rease in ICM densit y, such as that proposed for 
NGC4522 l| Kenney et al.l[2004l 'l. might have led to the dis¬ 
placement of the Hi and CO (J = 2 —^ 1) with respect to 
the stellar maxima without sufficient time for gas deficiencies 
or an Hi tail to develop. Such an increase in ram pressure 
might arise from shocks or ICM substructure which are more 
likely in a merging cluster, although the X-ray (XMM) data 
does not show evidence of an ICM density enhancement at 
or close to CGCG 97-079. For the ram pressure stripping 
archetype NGG4522 there is no indication that the highest 


density CO (J = 2 —>■ 1) was displaced from its currently 
observed location (in a disk surrounding the optical centre) 
during the proposed spike in ram pressure. Although there 
is kinem atic evidence that s ome Hi has fallen back subse¬ 
quently ll Kenney et al.l[2004l l. The argument that a similar 
order of magnitude spike in ram pressure is the cause the 
observed separation of the highest density CO (J = 2 —^ 1), 
HI and stellar components in CGCG 97-079 would therefore 
rely on this being the consequence of the shallower potential 
in CGGG 97-079 compared to NGC 4522. 


5 CONCLUDING REMARKS 

CGCG 97-079 is a large gas rich dwarf with a modest H l de¬ 
ficiency and a significant H 2 excess which displays complex 
morphology (and gas kinematics) at each of the wavelengths 
considered above. Taken together these clearly indicate the 
galaxy has suffered a strong perturbation. This perturbation 
probably coincided with a burst of star formation which evo¬ 
lutionary synthesis models indicate occurred ~ 10® yr ago. 

CGCG 97-073 and CGCG 97-079 display impressive 
70 and 100 kpc H q tails respectively in the deep 
iBoselli fc Gavaz^ (l2014h Ha image which mee t each other in 
projec tion NW of CGCG 97-079. Modelling in lCavazzi et al.l 
(I 2 OOII) suggests a high-speed tidal interaction, at the ver¬ 
tex of the Ha tails, subsequently enhanced the efficiency 
of ram pressure stripping of H l from both galaxies. The 
Ha tails of both g alaxie s only become visible in the deep 
iBoselli fc Gavaz^ |[20l3) Ha image from the tails’ vertex 
onwards. Prior to the high-speed tidal interaction, ram pres¬ 
sure stripping was insufficient to produce detectable Ha 
tails, i.e., both tails are detectable because the high-speed 
interaction “ loosened ” gas in the disks, which follow¬ 
ing stripping and ionisation now present as the Ha tails. 
The time scale for the high-speed interaction assuming 
a Vrei of 1000 kms~^ (which implies a current Pram = 
10 “^^ dyncm“^) is 10® yr in good agreement with the time 
sinc e the last burs t of sta r formation and the timescale from 
the ICavazzi et al.l (l200lll modeling. Based on the analysis 
in section im we concluded that H l, CO (J = 2 ^ 1) and 
stellar intensity maxima offsets (maxima offsets) are not the 
result of a high impact low velocity interaction with a near 
neighbour. We also investigated whether the maxima offsets 
could be the result a recent minor merger but did not find 
convincing evidence for this. A further possibility is the ram 
pressure stripping of gas has caused a shift in the cusp of the 
dark matter halo away from centre of the dark matter (DM) 
halo via drag forces. This effect was modelled for a medium 
mass (total mass = 10 ^° Mp,) dwarf sub ject to face-on ram 
pressure stripping by (ISmith et al.ll20l j l. Their model with 
a similar /gas (0.5) and ram pressure to that calculated for 
CGCG 97-079 predicts a shift in the cusp and stellar disk 
centres from the centre of the DM halo of < 0.5 kpc on time 
scales of 10® yr. But more critically in their models the gas 
within the truncation radius remains bound to the deep¬ 
est potential of the DM, i.e., this mechanism is unlikely to 
explain the high density offsets. The most likely remaining 
possibilities are that the offsets in neural and molecular ISM 
are the result of ram pressure and/or were produced by the 
high-speed interaction with CGCG 97-073. As we note in the 
introduction models and observations of spirals more mas- 

























11 



40.0" 


30.0" 


■B 20 . 0 " 


+ 20 ° 00 ' 00 . 0 " 


14.50s 


14.00s 


13.50s 13.00s 

RA (J2000) 


12.50s 


1.50 


1.35 


1.20 

1/1 

o 

. g 

1.05 M E 

c 

C 


0.90 


0.75 


0.60 


12.00s llh43mll.50s 






:U 




Figure 7. CGCG 97-079 SSDS g — i band image with the same CO { J = 2 —>■ 1) contours (blue) and Ho knot identifications as in 
Figure 2. 


sive than CGCG 97-079 indicate the estimated ram pressure 
of Pram = 10“^^dyncm“^ alone would not produce such 
offsets. Unlike ram pressure stripping, tidal interactions are 
known to be capable of affecting high density H l and molec¬ 
ular gas loc ated deep within the galactic gravitat ional po¬ 
tential well iIDuc fc Mirabellll994l : llono et ah 2005h and the 
relax ation time scale for Hi is > lx 10® yr I Holwerda et al.l 
[2 oT 3). High velocity "harassment” type fly-by interactions 
in higher mass spirals than CGCG 97-079 are expected to 
have only a minimal im pact on the old stellar component 
iIDuc fc Bournaudll2008l ). but it is unclear whether this also 
the case for lower mass spirals such as CGCG 97-079. For 
the CGCG 97-079 / CGCG 97-073 we do not know whether 
the interaction between them was a fly-by or the disks pen¬ 
etrated each other. 


In summary without clear evidence of the impact of the 
interaction on the old stellar disk (section |4.1 II and the un¬ 
certainly about whether the shallower gravitational poten¬ 
tial in CGCG 97-079, compared to the best studied cases, 
could allow ram pressure or a high-speed tidal interaction 
to produce the observed maxima offsets, we are unable to 
determine whether ram pressure or a high-speed tidal in¬ 
teraction was the principal cause of the observed maxima 
offsets. However ram pressure stripping is likely to be play¬ 
ing a significant role in the perturbation of lower density 
gas. Resolved ram pressure stripping modeling of the ISM 
in spirals with masses similar to CGCG 97-079 would greatly 
assist in resolving this question. 
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